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Pulsed laser excitation of a dense ultracold Cs vapor has been used to study the pairwise inter¬ 
actions between Cs atoms excited to np 3/2 Rydberg states of principal quantum numbers in the 
range n = 22 — 36. Molecular resonances were observed that correspond to excitation of Rydberg- 
atom-pair states correlated not only to the np 3/2 + np 3/2 dissociation asymptotes, but also to 
nsi /2 + (n + l)si/ 2 , nsi /2 -b n'fj, and (n — 4)fj + (n — 3)fj {j = 5/2, 7/2) dissociation asymptotes. 
These pair resonances are interpreted as arising from dipole-dipole, and higher-order long-range- 
interaction terms between the Rydberg atoms on the basis of i) their spectral positions, ii) their 
response to static and pulsed electric fields, and iii) millimeter-wave spectra between pair states 
correlated to different pair-dissociation asymptotes. The Rydberg-atom-pair states were found to 
spontaneously decay by Penning ionization and the dynamics of the ionization process were in¬ 
vestigated during the first 15 /rs following initial photoexcitation. To interpret the experimental 
observations, a potential model was derived that is based on the numerical determination of the 
eigenvalues and eigenfunctions of the long-range interaction Hamiltonian. With this potential model, 
which does not include adjustable parameters, all experimental observations could be accounted for, 
and the results demonstrate that long-range-interaction models provide a global and accurate de¬ 
scription of interactions in ultracold Rydberg gases and that they correctly account for, and enable 
the analysis of, phenomena as diverse as the formation of Rydberg macrodimers. Penning ionization 
in dense Rydberg gases, and Rydberg-excitation blockade effects. 


I. INTRODUCTION 

Long-range interactions between atoms in high Ryd¬ 
berg states are the origin of many scientific applications 
involving Rydberg excitation of ultracold atoms. The 
long-range van der Waals and resonant dipole-dipole in¬ 
teractions between high-n Rydberg states lead to effects 
such as the Rydberg excitation blockade mm that might 
be utilized in realizing scalable quantum gates for quan¬ 
tum information processing HHS], the evolution of ultra¬ 
cold Rydberg gases to ultracold plasmas [5] , as well as 
strongly nonlinear optical processes mm, e.g., to gen¬ 
erate single photons. While many of the recent experi¬ 
ments with ultracold Rydberg atoms rely on the inter¬ 
actions between the excited atoms, few experiments di¬ 
rectly probe the potential-energy structure of the inter¬ 
acting Rydberg gas. 

Rydberg-Rydberg interactions were first observed as 
density-dependent line broadening in the excitation of 
Cs atoms in a beam [T^] . The effects of these interac¬ 
tions became easier to observe in experiments on ultra¬ 
cold atoms because these enable much longer interaction 
times under conditions where the interatomic potential 
energy dominates the dynamics of the system. Asymmet¬ 
ric line shapes and interaction-induced shifts have been 
observed, e.g., in millimeter-wave spectra of transitions 
between Rydberg states m- Direct measurements of the 
long-range van der Waals and the resonant dipole-dipole 
interaction between two localized Rydberg atoms as a 
function of the interatomic separation were only reported 
recently oiin]. 
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We present experiments in which Cs Rydberg-atom 
pairs are detected as sharp resonances in laser-excitation 
spectra of ultracold, dense samples of Cs. The reso¬ 
nances are assigned and quantitatively analyzed using 
a model describing the interactions between Cs Rydberg 
atoms using long-range molecular potential curves that 
do not involve any adjustable parameters. Our experi¬ 
ments involve the pulsed excitation of pairs of interacting 
Cs Rydberg atoms with intense laser fields and enable 
the observation of molecular states connected to differ¬ 
ent pair-dissociation asymptotes. Whereas in our exper¬ 
iments only np Rydberg states are accessible by selec¬ 
tion rules for isolated atoms, Rydberg states can be ob¬ 
served as resonances near nsi/ 2 (fi +l)si/ 2 , ?t-Si/ 2 (?t- —3)fj, 
and (n — 4)fj(n — 3)fj asymptotes as a result of dipole- 
dipole, dipole-quadrupole m, and higher-order terms of 
the multipole expansion series. Our observations are re¬ 
lated to earlier work describing similar phenomena in 
ultracold Rydberg gases of rubidium HIHIO] and ce¬ 
sium mm at values of the principal quantum number 
n around 70. Compared to these previous investigations 
of molecular resonances, we study these resonances, re¬ 
ferred to as pair resonances hereafter, at lower principal 
quantum numbers in the range 22 < n < 36. Addition¬ 
ally to pair resonances observed in the immediate vicinity 
of dissociation thresholds, we also observe sharp spectral 
structures originating from local maxima in the molecu¬ 
lar potentials and from avoided crossings between differ¬ 
ent potential curves for n < 27. The positions, widths 
and relative intensities of these resonances are sensitive 
to the details of the molecular potentials at long range. 
We exploit this sensitivity to validate our Rydberg-atom- 
pair interaction model, which we then use to describe 
the Penning ionization processes that take place in our 
dense Rydberg gas on timescales of several hundreds of 
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nanoseconds, and the asymmetric broadening observed 
in millimeter-wave transitions between pairs of Rydberg 
atoms. 

This paper is organized as follows: The experimental 
apparatus and the measurement techniques are described 
in Section 11. The experimental results are presented in 
Section III. The Rydberg-pair-interaction model we de¬ 
veloped to account for the experimental observations is 
presented in Section IV. The main aspects of the paper 
are summarized in Section V and details concerning the 
calculations are presented in the appendix. 


II. EXPERIMENTAL 

All experiments are performed on ultracold samples of 
cesium atoms at a density of cm“^, a particle num¬ 

ber of ^10®, and a translational temperature of ~40/rK. 
The samples are released from a far-detuned crossed op¬ 
tical dipole trap (A = 1064 nm, P = 10 Watt, ~80^m 
1 /e^ radius) which is loaded from a magneto-optical trap 
(MOT) [23]. The transfer from the MOT to the dipole 
trap proceeds as follows: After a MOT-loading phase 
of ~80ms, the fiber laser used for the dipole trap is 
turned on and, simultaneously, the magnetic-field gra¬ 
dient is ramped up, compressing the MOT for 10 ms. 
The magnetic field is then turned off and the atoms 
are cooled in an optical-molasses phase for 5 ms. The 
MOT lasers are then turned off and the atoms are opti¬ 
cally pumped into the lower hyperfine component of the 
6 si /2 electronic ground state. The atoms are held in the 
dipole trap for 10 ms, during which time the noncap- 
tured atoms can leave the photoexcitation region. The 
power of the optical-dipole-trap laser is reduced to its 
minimum value of 300 mW after the cesium atoms in the 
optical dipole trap have been selectively pumped opti¬ 
cally into the upper hyperfine state of the 6 si /2 ground 
state, from where they are photoexcited in single-photon 
transitions to np 3/2 Rydberg states. The experimental 
cycle is repeated at a rate of 10 Hz. 

The excitation laser is a frequency-doubled pulse- 
amplified ring dye laser (Coherent 899-21) delivering 4.4- 
ns-long Fourier-transform-limited UV pulses (^319 nm, 
140 MHz bandwidth) with pulse energies of up to 100 /iJ. 
The laser is focused to a spot size of ^150/rm (1/e^ ra¬ 
dius) resulting in peak irradiances of up to lOOMW/cm^. 
The laser frequency is calibrated with a wavemeter (High 
Finesse WS6-200). We have also employed single-color 
two-photon excitation of nSi /2 ^ 6 si /2 transitions us¬ 
ing the pulse-amplified fundamental output of the ring 
laser at a wavelength of ^639 nm and a pulse energy of 
^0.5 mJ. 

Transitions are detected by monitoring the Cs“'" ion 
signal using a microchannel-plate (MCP) detector after 
state-selective pulsed-field ionization (PFI). For PFI, we 
employ a pulsed electric field with a rise time of 1 fis 
and maximal field strength = 1.25kV/cm. Only 

np Rydberg states with n > 26 can be field ionized 


under these conditions. We also observe the formation 
of ions on molecular resonances for Rydberg states with 
n < 26, but only after introducing some delay between 
laser and electric-field pulse. After PFI and before the 
next MOT-loading phase, the atom cloud is character¬ 
ized by saturated-absorption imaging. From the images, 
values for the particle numbers and densities are derived. 
On the atomic resonances the trap is fully depleted and 
the resulting ion signals strongly saturate our detection 
system, which is designed to be sensitive to single ions. 
The spectral positions of the atomic transitions are ob¬ 
tained from trap-loss spectra based on the saturated- 
absorption images. 

Millimeter-wave transitions between Rydberg states 
excited on atomic or pair resonances are measured using 
a millimeter-wave source consisting of a swept frequency 
synthesizer (Wiltron, Model 6769B) connected to an 
amplifier-multiplier chain (Virginia Diodes, WR6.5AMC, 
tunable range 110-170GHz). Typical interaction times 
with the millimeter-wave radiation are 1 fis and the tran¬ 
sitions are detected using state-selective PFI and sep¬ 
arating initial and final Rydberg states in the Cs^-ion 
time-of-flight (TOF) traces. 

Regular measurement and compensation of the stray 
electric and magnetic fields during the data acquisition 
period ensured that the magnitude of the stray elec¬ 
tric and magnetic fields never exceeded 20 mV/cm and 
20 mG, respectively, which is sufficiently low that these 
fields do not influence the behavior of the Rydberg states 
we study. The electric field was compensated by measur¬ 
ing and minimizing the quadratic Stark effect of high 
^P 3/2 Rydberg states [23| and the magnetic field was 
compensated by measuring the Zeeman splitting of the 
F' = A F = 3 clock transition of the 6 si /2 state using 
microwave spectroscopy. 


III. EXPERIMENTAL RESULTS AND 
DISCUSSION 

A. Excitation of dipole-coupled pair states 

The experimental spectra presented in this section 
were recorded by monitoring the Gs“'"-ion yield resulting 
from the spontaneous ionization of the Rydberg atoms 
taking place in the first 5 /iS following photoexcitation in 
the dense ultracold sample as a function of the UV-laser 
frequency. Because of a Gooper minimum in the pho¬ 
toexcitation cross section from the 6 si /2 ground state of 
Cs to the npi /2 (ePi/ 2 ) states (continuum) with center 
located just above the ionization threshold [ 21 |, transi¬ 
tions to npi /2 Rydberg states are more than 500 times 
weaker than to the corresponding np 3/2 Rydberg states. 
Consequently, the experiments focus on the np 3/2 states. 
To facilitate the comparison between spectra recorded in 
the vicinity of different np 3/2 6 si /2 transitions, the 

origins of the frequency scale of the different spectra pre¬ 
sented in the following are always set at the positions of 
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the atomic resonances. 

Fig .0 compares two spectra of the 31p3/2 ^ 6 si /2 
transition recorded at a ground-state-atom density of 
^ 10 ® atoms/cm^ (dashed line, sample obtained by re¬ 
leasing the atoms from the MOT and letting them ex¬ 
pand) and of ~5 • 10^^ atoms/cm® (full line, sample ob¬ 
tained after release from the optical dipole trap). In 
both spectra, the atomic 31p3/2 ^ 6 si /2 transition is 
strongly saturated. The spectrum measured at lower 
atom density has been scaled up to match the wings of 
the 31p3/2 6 si /2 transition recorded at higher atom 

density. In the latter spectrum, an additional resonance 
is observed at a detuning of —1200 MHz. The spectral 
position of this resonance cannot be assigned to any 
isolated-atom transition in Cs but matches closely the 
asymptotic energy of the 31 si/ 232 si /2 pair state if one 
assumes a two-photon excitation, for which the actual 
detuning is 2 x (—1200MHz). The resonance has a red- 
degraded asymmetric line shape and a full width at half 
maximum (FWHM) of 2 x 250 MHz. 

In the absence of interactions between the Rydberg 
atoms, the single-color two-photon excitation of the 
3 Isi/ 232 si /2 ^ 6 S 1 / 26 S 1/2 transition is prohibited be¬ 
cause pure s' •<— s transitions are forbidden by the La- 
porte rule. However, under our experimental conditions, 
long-range interactions between the Rydberg atoms can 
lead to state mixing, especially between energetically 
close-lying dipole-dipole-coupled pair states such as the 
nsi/ 2 {n + l)si /2 and the np 3 / 2 ?T-P 3/2 states, as schemati¬ 
cally depicted in Fig. The molecular shifts AEi^t in¬ 
duced by the long-range interactions are expected to be 
small and resonances occur close to the nsi/ 2 {n + l)si /2 
asymptotes. The n dependence of the pair-state asymp¬ 
tote positions with respect to the atomic np 3^2 resonance 
positions can be calculated as 

A / _ AE{n) _ ~ ^riSi /2 ~ ^{n+l)si /2 

with 

_ hRcs . , 

In Equation i?Cs = 3.289828299(20) • 10® MHz is the 
Rydberg constant of Cs and Si^ (n) are the I- and weakly 
j- and n-dependent quantum defects. The n dependence 
of the quantum defects can be expressed with the Ritz 
formula [25] 


Si^ {n) = 61 ^ ( 00 ) 


(n- i5q(oo))2' 


( 3 ) 


The parameters Si. { 00 ) and ai. are taken from the liter¬ 
ature [55]. The detuning A^(n) of the pair asymptotes 
decreases rapidly in the range from n = 20 to n = 42, 
where the detuning changes sign. In the range of n val¬ 
ues investigated, Ai/(n) approximately scales as (see 
fit shown as blue line in Fig.[^). Because the frequency 
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FIG. 1. UV-laser spectra near the 31p3/2 ■«— 6 si /2 transition 
of Cs at ground-state-atom densities of around 10® atoms/cm® 
(dashed trace) and of 5 • 10®® atoms/cm® (solid trace). 



FIG. 2. (Color online) a) Simplified and schematic dia¬ 
gram of the interaction between the nsi/ 2 (n + l)si /2 and the 
’^P 3 / 2 np 3/2 Rydberg-atom-pair states. The blue thick and 
the thin, dashed lines indicate strong and weak transitions, 
respectively, from the 6 si /2 ground state. Two-photon tran¬ 
sitions are indicated by two parallel arrows, b) Energetic 
separation between the two pair states as a function of the 
principal quantum number n. The blue curve corresponds to 
the function an”®'®® obtained by fitting the exponent to the 
energies of the dissociation asymptotes of nsi/ 2 (»^+l)si /2 pair 
states relative to the np 3 / 2 U.p 3/2 asymptotes. 


required for the two-photon excitation of the pair reso¬ 
nance is detuned from the atomic np 3/2 t— 6 si /2 reso¬ 
nance, the two-photon resonance can only be observed 
at the high laser power achievable with a pulse-amplified 
laser system. It is also necessary to have a high ground- 
state-atom density to have a sufficient number of atom 
pairs at distances where the dipole-dipole interaction is 
strong enough to induce the state mixing illustrated in a 
simplified 2-level model in Fig. i®)- No pair resonances 
could be observed in spectra recorded with a 100 mW 
continuous-wave UV laser under otherwise similar exper¬ 
imental conditions m- 

Rydberg-excitation spectra in the vicinity of np 3/2 t— 
6 si /2 transitions for n = 22 — 36 are displayed in Fig.j^ 
At each n value, we observe resonances on the low- 
frequency side of the atomic transitions that can be at¬ 
tributed to the two-photon excitation of two atoms to 
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interacting pair states correlated to the nsif 2 {n + l)si /2 
asymptotes. The positions of these asymptotes calcu¬ 
lated with Eq. (|^ are marked by vertical gray lines in 

Fig.[| 

For all n values for which measurements could be per¬ 
formed, i.e., for n = 22 — 36, pair resonances are observed 
within 2 X 75 MHz of the positions of the riSi/ 2 {n + l)si /2 
asymptotes. The deviations are displayed on a magnified 
scale in the inset of Fig.|^ (open circles). In most cases, 
the deviations are negative, which indicates that the pair 
resonances are typically located on the low-frequency side 
of the asymptotes, which is characteristic of attractive 
molecular potentials. The measurement at n = 25 repre¬ 
sents an exception and reveals a pair resonance located 
2 X 75 MHz above the 25si/226si/2 asymptote. These 
shifts from the asymptotes, though small, are real and 
will be discussed in more detail in Section IV. 

The increasing detuning of the pair resonances from 
the atomic transitions (Fig.[^ and the decreasing inter¬ 
action strength with decreasing n values lead to signifi¬ 
cantly smaller pair-state signals at lower n values. This 
effect is compensated by using higher pulse energies in 
the excitation, as indicated on the left-hand side of Fig. 

At the lowest n values for which the Rydberg-atom-pair 
resonances could be observed {n = 22 — 26 in Fig. [^, we 
also observe additional structures on their high-frequency 
side. These structures are discussed in more detail in Sec¬ 
tion IV. The dipole-quadrupole pair resonances presented 
in Ref. [TB] are located at smaller detunings and are off 
the vertical scale of the panels of Fig.[^ 

The close agreement between the positions of the 
nsi/ 2 (?T-+l)si /2 asymptotes and the molecular resonances 
near the np 3/2 6 si /2 atomic transitions in high-density 

samples (see Figs.[^ and strongly suggests that the 
mechanism responsible for their observation is the dipole- 
dipole-coupling mechanism illustrated schematically in 
Fig.[2}i). Further evidence for this mechanism was ob¬ 
tained from the measurement of the Stark effect and of 
millimeter-wave transitions discussed in the next subsec¬ 
tions. 


B. Stark effect of the nsi/2(u-|- l)si/2 and UP3/2UP3/2 
pair states 

The atomic s and p Rydberg states of Cs exhibit 
quadratic Stark shifts in small external electric fields be¬ 
cause these states have large quantum defects and are 
well separated from other Rydberg states of opposite 
parity. Therefore, when a small external electric field 
of strength F is applied, one expects the positions of 
the np 3/2 6 si /2 transitions (and hence the positions of 

the np 3 / 2 ?T-P 3/2 asymptotes) to be shifted by 
and the positions of the nsi/ 2 (n-l-l)si /2 asymptotes to be 
shifted by \F'^ [\as^j^{n) + 1)] in UV spec¬ 

tra of the type depicted in Figs.[^and[^ The polarizabil¬ 
ities ai- are negative and much larger for P 3/2 Rydberg 
states than for Si /2 Rydberg states of the same n value. 
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FIG. 3. (Color online) Pulsed-UV-laser spectra of transitions 
from the 6 si /2 ground state of Cs to np 3/2 Rydberg states 
with n between 22 and 36 (solid black traces). The Cs''"-ion 
signal is plotted as a function of the detuning from the atomic 
^P 3/2 6 si /2 transition. The spectral positions correspond¬ 

ing to the asymptotic energy of nsi/ 2 (u + l)si /2 pair states 
calculated from experimental quantum defects are indicated 
by grey vertical lines. Simulated spectra are shown as red 
dashed lines. The pulse energies used to record the different 
experimental spectra and the relative intensity scaling factors 
of the experimental spectra are given on the left-hand side of 
each panel. 


Consequently, a measurement of the Stark shifts of the 
pair resonances at low fields, for which the Stark shift is 
quadratic, enables one to distinguish between pair states 
of rip 3 / 2 np 3/2 and nsi/ 2 ('u + l)si /2 character. The dif¬ 
ferent polarizabilities of these two pair states can also be 
used to bring them into resonance with electric fields and 
to observe effects of resonant dipole-dipole interactions, 
which scale as jF? O UHl US]- For the measurements 
relevant for this investigation, the dipole-dipole interac¬ 
tion between np 3 / 2 np 3/2 and nsi/ 2 (u-l- l)si /2 pair states 
is nonresonant and of shorter range (scaling with 1 /i?®) 
so that the measurement of the pair-state Stark shifts ne- 
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FIG. 4 . Spectral positions of the nsi/2(n + l)si/2 pair res¬ 
onances relative to the atomic np3y2 states for values of the 
principal quantum number n in the range between 22 and 36 . 
The open circles represent the peak positions extracted from 
measured spectra and the black line corresponds to the rela¬ 
tive positions of the asymptotes calculated from experimental 
values of the quantum defects of Cs using Eq. Q. The exper¬ 
imental (simulated) interaction-induced shifts AEint relative 
to the nsi/2(n-|- l)si/2 pair-dissociation asymptotes are given 
in the inset as open circles (dashed line). 


cessitates high laser intensities and high atom densities. 

We have measured Rydberg excitation spectra of Cs in 
the vicinity of the 32p3/2 t— 6 si /2 transition in weak ex¬ 
ternal electric fields between OV/cm and 12V/cm. The 
electric-field strength was varied by applying different 
electric potentials to sets of electrodes surrounding the 
ultracold atoms, as described in Ref. [53]. The measured 
spectra have the same general appearance as the spectra 
depicted in Fig.|^ but the strong atomic and weak pair 
resonances are shifted by the quadratic Stark effect, and 
the atomic resonance splits into the two magnetic compo¬ 
nents with \mj\ = 1/2 and 3/2. All measured positions 
with respect to the field-free 32p3/2 ^ 6 si /2 atomic tran¬ 
sition frequency are plotted as a function of the applied 
electric field in Fig.[^ where the black and red dots cor¬ 
respond to the atomic and pair resonances, respectively. 
The observed quadratic Stark shift and splitting of the 
32p3/2 Rydberg state was used to calibrate the applied 
electric field by comparison to the calculated Stark shifts 
of the \mj\ = 3/2 (dashed line) and \mj\ = 1/2 (dotted 
line) components of the 32p3/2 state. The Stark shift of 
the 32 si/ 233 si /2 pah state is much smaller and matches 
the average of the calculated Stark shifts of two non¬ 
interacting atoms in the 32si/2 and 33 si /2 states (red 
line in Fig.j^ within the experimental uncertainty. The 
32 p 3 / 232 p 3/2 and the 32 si/ 233 si /2 asymptotic energies 
are degenerate at an electric field of « 7 V/cm. At higher 
electric-field strengths, we observe the 32 si/ 233 si /2 reso¬ 
nance on the high-frequency side of the transitions to the 
32p3/2|mj| = 1/2,3/2 states. The agreement between 
measured and calculated Stark shifts indicates that the 
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FIG. 5 . (Golor online) Stark shifts of the atomic 32 p 3/2 
6si/ 2 transition of Gs (black points) and of the 32 si /2 33si/2 t— 
6S1/26S1/2 two-photon pair resonance (red points). The calcu¬ 
lated quadratic Stark shifts of the atomic and the pair (aver¬ 
aged atomic Stark shifts of the two separate Rydberg atoms) 
resonances are given as dashed and dotted black and solid red 
lines, respectively. 


Rydberg atoms excited at the positions of the pair res¬ 
onances have predominantly S 1/2 character and that the 
pairs are excited close to their asymptote, at long range. 


C. Millimeter-wave spectroscopy of transitions 
between Rydberg-atom—pair states 

The ability to observe well-resolved transitions to 
nsi/ 2 (n + l)si /2 dipole-dipole-coupled pair states with 
good signal-to-noise ratio offers the possibility to further 
characterize these states by studying transitions to neigh¬ 
boring dipole-dipole-coupled states using millimeter- 
wave radiation. In a similar fashion microwave transi¬ 
tions between pairs of Rb Rydberg atoms were studied 
by Y. Yu et al. |30j . 

At large internuclear distances, one expects the most 
intense transitions to be those to pair states located near 
nsi/ 2 nipj and n 2 Vjin+l)si /2 asymptotes, where the val¬ 
ues of the principal quantum numbers ni and 712 are 
close to the value of n. Our investigation focused on 
the 32 si/233si/ 2 pair state and on the millimeter-wave 
transitions to the pair states near the 32si/233p3/2 and 
31p3/233si/2 asymptotes, which lie 129611MHz above, 
and 157962 MHz below the 32si/233si/2 asymptote, re¬ 
spectively (see Fig.|^)). The detection of the millimeter- 
wave transitions was achieved by delayed pulsed-field ion¬ 
ization with a slowly rising electric-field ramp (rise time 
1/rs) and monitoring the Cs+-ion time-of-flight (TOF) 
spectra. To calibrate the TOF spectra, experiments were 
first performed in low-density samples following single¬ 
photon laser excitation to the np 3/2 {n = 31,33) and 
nonresonant single-color two-photon excitation to the 
n'si /2 {n' = 32,33) Rydberg states of Cs. The result- 
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FIG. 6. (Color online) a) Schematic energy diagram of the 
32 si/233si/ 2 ■<— 6S1/26S1/2 pair state excitation (left-hand 
side, blue arrows) and the single-color two-photon excitation 
of nsi/2 ■‘r- 6si/ 2 transitions (right-hand side, red arrows). 
Millimeter-wave transitions are indicated with black arrows 
in the energy-level diagrams. Experimental millimeter-wave 
spectra of the 33p3/2 33si/2 (b) and 32 si /2 —>■ 31 p 3/2 (c) 

transitions recorded following excitation of the 32si/233si/2 
pair state (solid black traces) and following excitation of the 
atomic nsi/2 6si/2 transitions (dashed black traces), re¬ 
spectively. The spectra represent the ratio of ions detected in 
the integration window of the relevant final Rydberg states 
over the ions detected in the integration windows correspond¬ 
ing to both the final and initial Rydberg states (see Fig.[^. 
The red dashed lines are spectra simulated with a model tak¬ 
ing into account all relevant Rydberg-Rydberg-pair poten¬ 
tials and considering the Penning-ionization dynamics, as ex¬ 
plained in more detail in the text (see Section IV D). 


ing Cs’*' TOF traces are displayed in the lower panel of 
Fig.0 and reveal distinct maxima for the four Cs Ryd¬ 
berg states. This state selectivity makes it possible to 
record millimeter-wave spectra of the 33p3/2 ^ 33si/2 
and 32 si/ 2 31p3/2 transitions by monitoring the Cs+- 

ion signal at the corresponding TOF positions as illus¬ 
trated by the spectra displayed as dotted lines in Fig. 
and c), respectively. The lines in these spectra have a 
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FIG. 7 . (Color online) Cs"'"-ion time-of-flight (TOF) distribu¬ 
tions measured using a pulsed extraction field of 1 . 25 kV/cm 
with a 1 /rs rise time. Top panel: TOF distributions following 
the two-photon excitation of the 32si/233si/2 pair state with¬ 
out (red trace) and with millimeter-wave radiation (mmW) 
with the frequency tuned to the 33p3/2 33si/2 transition 

(brown trace). Bottom panel: TOF distributions following 
the one-photon excitation of the atomic 33p3/2 (black trace) 
and 32 p 3/2 (blue trace) states, and the two-photon excitation 
of the atomic 33si/2 (violet trace) and 32 si /2 (turquois trace) 
states. The nsi/2 states have been excited from the 6si/2 
ground state in single-color two-photon transitions. Hatched 
areas indicate the integration windows used to extract the 
spectra shown in Fig.|^). 


near-Fourier-transform-limited width of 2 MHz. 

The Cs^-ion TOF spectrum obtained at the position 
of the 32si/233si/2 pair resonance and a delay of 1 /rs be¬ 
tween the laser and the electric-field pulse is displayed as 
red trace in the upper panel of Fig.[^ It consists of an in¬ 
tense sharp peak at early times which is attributed to Cs'*' 
ions produced before the electric-field pulse (the dynam¬ 
ics of the ionization process is discussed in Section HID), 
and a broad distribution with only weak maxima at the 
positions expected for the pulsed-field ionization of the 
atomic 32si/2 and 33si/2 states. 

Compared to the TOF spectra recorded for the non¬ 
interacting Cs Rydberg gas, the state selectivity of the 
pulsed-field-ionization process is reduced. The reduc¬ 
tion is attributed at least in part to adiabatic traver¬ 
sals of the curve crossings between different pair states, 
e.g. between the 32si/233si/2 and 32p3/232p3/2 states 
(see Fig.[^ before and as the electric field is ramped up. 
A similar observation was made by Han and Callagher 
in their study of the field-ionization dynamics in an ul- 
tracold Rb Rydberg gas (see in particular their Fig. 4 in 
Ref. [U). The reduced state selectivity of the field ion¬ 
ization causes a millimeter-wave-radiation-independent 
background signal at the TOF positions correspond¬ 
ing to the 31p3/2 and 33p3/2 Rydberg states, which 
makes the recording of millimeter-wave spectra of the 
(31p3/233si/2j 32si/233p3/2) 32 si/233si/2 transitions 

difficult. However, comparing the Cs^-ion TOF spec¬ 
trum recorded with and without millimeter-wave radia¬ 
tion under otherwise identical conditions reveals a weak 
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enhancement of the 33p3/2 PFI signal and a weak deple¬ 
tion of the 33 si/ 2 PFI signal when the millimeter-wave 
frequency is tuned to the 33p3/2 ^ 33 si/ 2 resonance (see 

Fig.0. 

The ion-signal ratios at the TOF positions of the 
33p3/2 and 33si/2 states and at the positions of the 31p3/2 
and 32 si/ 2 states were used to record millimeter-wave 
spectra of the transitions from the 32 si/233si/ 2 pair state 
to regions near the 32si/233p3/2 and 31p3/233si/2 asymp¬ 
totes, respectively. These spectra are presented as solid 
lines in Fig.|^ The observation of a millimeter-wave- 
dependent PFI signal corresponding to the 31p3/2 and 
the 33p3/2 Rydberg states represents a further experi¬ 
mental verification of the 32si/233si/2 nature of the laser- 
excited pair resonance. Compared to the millimeter-wave 
spectra recorded in the noninteracting Rydberg gas, the 
transitions are broadened and degraded toward lower 
(higher) frequencies in the case of the 32si/233si/2 —> 
31p3/233si/2 (32s3y233p3/2 ^— 32 S 3 / 233 S 3 / 2 ) transition. 
Although the large background signal and the poor 
signal-to-noise ratio of the millimeter-wave spectra of the 
pair states prevents a quantitative analysis, the spectra 
enable us to draw the following qualitative conclusions: i) 
The weak line broadenings and frequency shifts compared 
to the transitions recorded in the noninteracting Ryd¬ 
berg gas suggest that the Rydberg-atom pairs contribut¬ 
ing to the spectra have a large interatomic separation, ii) 
The red (blue) shift of the transitions to the 31p3/233si/2 
(32si/233p3/2) pair state reveals that the potential curves 
of the n'p 3 / 2 nsi /2 pair states are less strongly attractive, 
at long range, than is the case for the 32si/233si/2 pair 
state, despite the fact that the dipole-dipole interaction 
is resonant, because nSi/ 2 ''T- 2 Pj and nipjnsi /2 are degen¬ 
erate when ni = n 2 . This aspect is discussed in more 
detail in Section IV D. 

The observation of a prompt ionization signal further 
shows that an important fraction of the initial population 
in the 32si/233si/2 pair state decays by ionization on the 
sub-microsecond timescale and therefore does not con¬ 
tribute to the millimeter-wave spectra. The millimeter- 
wave spectra can be adequately modeled (see red dashed 
traces in Fig.[^ by taking into account the motion of the 
Rydberg-atom pairs on their potential-energy curves, as 
explained in Section IV D. 


D. Penning ionization of the 32si/233si/2 pair state 

The spontaneous ionization of Rydberg atoms in the 
gas phase is the first step in the evolution of a dense 
Rydberg gas towards a plasma. Spontaneous ioniza¬ 
tion of a Rydberg gas was first reported by Vitrant et 
al. [ 35 ], who pointed out that it reduces the lifetimes of 
Rydberg states in dense samples. Interest in the ion¬ 
ization process was renewed with experiments on ultra¬ 
cold atoms which provided detailed information on the 
evolution from an ultracold Rydberg gas to an ultra¬ 
cold plasma [ 33 ] . This evolution can be described by 


an avalanche-ionization model [34j that involves an ini¬ 
tial seed ionization process, during which ions are pro¬ 
duced in the ultracold gas and the free electrons escape. 
Once the positive space charge of the ions is sufficiently 
high to trap low-energy electrons, the remaining Ryd¬ 
berg atoms are rapidly ionized by collisions with the 
electrons. For the seed-ionization process, experiments 
pointed at the importance of a small fraction of hot Ry¬ 
dberg atoms and blackbody-radiation-induced photoion¬ 
ization [33] . static orbital overlap between close-lying 
Rydberg-atom pairs, sequences of quasiresonant dipole- 
dipole transitions PI3S], and interaction-induced motion 
leading to Rydberg-Rydberg collisions and Penning ion¬ 
ization [36ll38] . 

The dynamics leading to the seed ionization depends 
on the density of the excited Rydberg atoms, on the Ry¬ 
dberg state and on the details of the interaction between 
the Rydberg atoms. At high n values (n « 100) and high 
Rydberg-atom densities (> 10®cm“^), the complete ion¬ 
ization of the Rydberg gas occurs in less than 100 ns, 
as reported in Ref. |31| and confirmed in the present 
work (not shown). The theoretical understanding of this 
fast ionization process is still incomplete, but a classi¬ 
cal simulation [39] sugggests that many-body interactions 
and excitation of closely-spaced pairs of Rydberg atoms 
play an important role under these conditions. At lower 
densities and n values, Penning-ionization processes [331 - 
iMiiiniiii] involving two Rydberg atoms, one atom un¬ 
dergoing a transition to a lower state and the other be¬ 
ing ionized, have been identified as the dominant seed- 
ionization mechanism. 

To characterize the processes leading to the strong 
prompt Cs“'" ionization signal in our studies of millimeter- 
wave transitions between Rydberg-pair states (see up¬ 
per panel of Fig.[^, the ionization dynamics of the 
32si/233si/2 pair state was studied at a time resolution of 
100 ns limited by the temporal resolution of our Cs^-ion- 
TOF spectra. The experiments consisted in monitoring 
the yield of Cs’*' ions produced spontaneously during the 
delay time between the laser excitation pulse and the ap¬ 
plication of the pulsed electric field as a function of the 
UV-laser frequency. The spontaneously formed ions are 
accelerated from the beginning of the electric-field ramp 
used for the field-ionization and thus arrive first on the 
detector, see the integration-windows labeled “prompt 
ions” in Fig. [^ 

The results of this investigation are presented in 
Fig.||r), which compares spectra of the 32si/233si/2 pair 
resonance recorded for several delay times between 200 ns 
and 9.9 /iS. At the shortest delay times, the peaks in the 
prompt-ion signal are observed « 2 x 300 MHz below the 
32si/233si/2 asymptote and are very weak. The intensity 
of the resonance grows with increasing delay time and its 
position gradually shifts toward the 32si/233si/2 asymp¬ 
tote, and stabilizes at « 2 x 75 MHz below this asymp¬ 
tote at delay times beyond 1 ^s. This deviation from the 
asymptotic pair energy is the interaction-induced shift 
AAint discussed in section jlll Aj and shown in the inset 





FIG. 8. (Color online) a) Spontaneous-ionization signal as a function of the detuning of the UV-laser frequency from the 
32 p 3/2 -^681/2 transition and of the delay between laser pulse and pulsed ion-extraction field (solid colored traces, values for 
the delay time given in the upper right inset). The spectral position corresponding to the asymptotic energy of the 32si/233si/2 
pair state is shown by a vertical grey bar. Simultaneously recorded spectra of the total ion signal including ions resulting 
from pulsed-field ionization are shown as dashed colored traces on the same scale. The interaction-induced shift (i.e., twice 
the detuning of the resonance peak from the asymptotic energy of the 32si/233si/2 pair state) is given in the upper left inset 
and the fractions of spontaneous ion to the total-ion signals at the detunings indicated by red, black, and blue lines on the 
horizontal axis are displayed in the lower left inset, as functions of the delay time. Experimental values are shown as dots 
with error bars and simulated quantities as solid lines, b) Spectra simulated for the same experimental parameters as in panel 
a), as described in Sec. IV C. The inset shows the potential-energy functions (the two = 0 (green, blue, and black) and the 
n = 1 (red) potentials associated with the 32si/233si/2 asymptote are shown in color) of pair states below the 32p3/232p3/2 
asymptote as a function of the internuclear separation. 


of Fig. 1^ The time dependence of the resonance position 
is displayed in the upper left inset of Fig.|^ 

Because the total density of Rydberg atoms in our 
experiments is still very low (on the order of 1 • 10® 
atoms/cm®), the experimental observations can be qual¬ 
itatively understood as resulting from the purely pair¬ 
wise Penning-ionization mechanism. Pairs of Rydberg 
atoms are excited at well-defined initial separations given 
by the detuning of the excitation laser from the dissoci¬ 
ation asymptote and are accelerated on the attractive 
potential-energy curves. The first pairs that ionize are 
those characterized by the strongest shifts, i.e., those that 
correspond to the shortest internuclear distances and are 
thus subject to the strongest attractive interactions. As 
the delay time increases, atom pairs subject to smaller 
shifts and thus weaker attractive interactions start con¬ 
tributing to the prompt-ion signal in the TOF spectra. 
These pairs have larger initial internuclear separations 
and require longer times before the motion on the at¬ 
tractive potential leads to ionizing collisions. A more 


quantitative analysis of this process is presented in Sec¬ 
tion IV C. 


E. Spectral structures originating from pair states 
on the high-frequency side of the np 3 / 2 np 3/2 
pair-dissociation asymptotes 

The interaction between the np 3 / 2 ?T-P 3/2 and the 
nsi/ 2 (n + l)si /2 pair states schematically represented 
in Fig. leads to predominantly attractive potential- 
energy functions correlated to the nsi/ 2 ('n -I- l)si /2 
dissociation asymptotes and repulsive potential-energy 
functions correlated to the np 3 / 2 np 3/2 asymptotes for 
the considered range of principal quantum numbers n. 
Whereas the former lead to isolated resonances on the 
low-frequency side of the strong atomic resonances (see, 
e.g., Figs.[^ an d§ , the latter lead to dense and com¬ 
plex spectral structures on the high-frequency side of 
the atomic resonances, which are partially resolved at 
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FIG. 9 . (Color online) Spectral structures observed on the high-frequency side of the np3y2n-P3/2 (n = 22 — 26 ) dissociation 
asymptotes in the UV excitation spectrum of Cs. The experimental spectra (full black lines, identical in all panel) represent the 
Cs+ -ion signal as a function of the detuning of the UV-laser frequency from the np3/2 ^ 6si/2 transition. The pulse energies 
used to record these spectra and the relative intensity-scaling factors are indicated below the experimental traces on the left- 
hand side of the left panel. The panels also display simulated spectra (dashed lines) obtained by considering all interaction 
terms in the multipole expansion up to dipole-dipole (left, gray), quadrupole-quadrupole (middle, red), and octupole-octupole 
(right, blue) interaction terms. 


the lowest n values studied (n = 22 — 26). The corre¬ 
sponding spectra are displayed on an expanded scale in 
Fig.i where they are compared to numerical simulations 
performed using the interaction model presented in Sec¬ 
tion IV A. 


The observed spectral patterns, which are perfectly re¬ 
producible in the experiments, are more complex than 
those observed on the low frequency side and do not re¬ 
veal any obvious scaling with the principal quantum num¬ 
ber. They originate from the interactions of the repulsive 
’^P 3 / 2 ’^P 3/2 pair states with attractive (n — 4)fj(n — 3)fj 
pair states which have their dissociation asymptotes lo¬ 
cated just above the np 3 / 2 np 3/2 asymptote. The nu¬ 
merical simulations of the spectra in these regions in¬ 
dicate that the spectral structures do not only arise 
from an accumulation of spectral intensity in regions 
where the potential-energy functions are flat as in case of 
the nsi/ 2 {n + l)si /2 resonances, but also from a reduc¬ 
tion of spectral intensity centered around avoided cross¬ 
ings between states correlated to the np 3 / 2 rrp 3/2 and 
(n — 4)fj(n — 3)fj dissociation asymptotes. This aspect 
will be discussed in more detail in Section IV B. 


IV. POTENTIAL MODEL FOR INTERACTING 
RYDBERG-ATOM PAIRS AND 
INTERPRETATION OF EXPERIMENTAL 
RESULTS 


A. Model of long-range potential-energy curves 


The potential model we have developed to analyze the 
results presented in Section III was briefly introduced in 
Ref. [H] and is presented in detail here. We describe the 
long-range interaction between the two neutral atoms A 
and B in the frame of the nonrotating molecule AB using 
a multipole expansion of the form [32] 


binter 

oo +^< 




^A/B —1 OJ — —L< 


QL,u>ir) = 


47r 

2L + 1 


r^YL,uiir) 


( 4 ) 

( 5 ) 


[La + Tb)! 


\/ {La + w)!(La ~ w)!(Lb + a;)!(LB ~ w)! 

( 6 ) 
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where fi is the position of the Rydberg electron bound 
to atom i with respect to its core, R the interatomic 
separation, QLiuifi) the atomic multipole of L^-th order 
(L< being the smaller of two Li values), and w is the 
quantum number associated with the projection of the 
angular momentum L< onto the molecular axis. This 
expansion is valid in the absence of retardation effects, 
which are negligible in the systems discussed here, and 
for interatomic separations exceeding the LeRoy-radius 
Rlr = 2 [IS], a criterion specifying 

the internuclear distance beyond which orbital overlap 
between the two atoms becomes negligible. The atomic 
multipoles QLiuii^i) are expressed in terms of spherical 
harmonics (see Eq. (1^). To account for the experi¬ 
mental observations, we found it sufficient to restrict the 
multipole expansion to contributions up to the octupole- 
octupole (La = 3, Lb = 3) term. The Hamiltonian 

H = (7) 

A B 

consists of the Hamiltonians Hq ' of the two isolated 
atoms taking into account the quantum defects of all 
states with I < 4 (from Refs. HHm]) and the truncated 
multipole expansion term Einjer- 

To obtain the adiabatic potential-energy curves of 
the interacting Rydberg-atom pairs at distances larger 
than the LeRoy radius, H is determined in matrix 
form using the product basis Iua^aJawAi = 

lyAjAWA, 7 BjBWB), following the approach described in 
Refs. [TBH5U] . Here, wa.b represent the projections of 
the Rydberg-electron angular momentum onto the in¬ 
ternuclear axis, with wa + wb = We neglect the 
hyperfine structure of the Rydberg states because the 
largest hyperfine splittings for the low-n Si /2 states are 
still more than one order of magnitude smaller than our 
laser bandwidth [23] . The angular and radial parts of the 
matrix elements of H can be separated and the angular 
matrix elements are calculated analytically, whereas the 
radial matrix elements are obtained by numerical inte¬ 
gration US) using experimental quantum defects [231 Slj 
and an ab-initio model potential [46] . The details of the 
calculation of the matrix elements are presented in the 
appendix. 

The quantum number LI associated with the projec¬ 
tion of the total electronic angular momentum onto 
the molecular axis is a good quantum number. As 
pointed out in Ref. m, although coupling terms with 
(_l)iA-i-LB _ dipole-quadrupole cou¬ 

pling term, do not conserve the parity of the electronic 
wave function, they nevertheless lead to observable in¬ 
teractions because the coupling of quasi-degenerate rota¬ 
tional levels ensures the conservation of the total (rovi- 
bronic) parity. In Ref. [TS], we erroneously also stated 
that such terms can be accompanied by a mixing of 
states of gerade (g) and ungerade (u) electronic sym¬ 
metry in homonuclear dimers. The g-u mixing in Cs 
is actually induced by hyperfine interactions, which are 
neglected in our model. The eigenvalues of the Hamil¬ 


tonian H are determined in an unsymmetrized basis in¬ 
dependently for the different values of H = (0,1,2,3) 
that are optically accessible, and for different values of 
the internuclear separation R. The molecular eigenfunc¬ 
tions l^'(R)) = Ec 7 aJa<.^a. 7 bJbc^b(^) l7AjAWA,7BiBWB) 
and the molecular eigenenergies E^},{R) are determined 
as a function of the interatomic separation R. Exem¬ 
plary potential-energy curves describing the long-range 
interaction between two Cs Rydberg atoms are displayed 
in Figs.[lO|and[TT| 

The density of Rydberg-atom pair states is very high 
because of the large number of pair dissociation asymp¬ 
totes and because several potential-energy functions cor¬ 
relate to each asymptote. Only a few of the multi¬ 
tude of molecular pair states are accessible optically in 
our experiments. The excitation probability p^{R) of a 
molecular eigenstate |4>) at a pair separation R is pro¬ 
portional to the square of the mixing coefficients of the 
np 3 / 2 np 3/2 basis states. The optically accessible states 
\(l)) (for np 3 / 2 np 3/2 we take n,l = l,jA = Jb = 3/2) are 
conveniently expressed in the laboratory-fixed frame as 

I/*) = ^ ajAmAOBlTlB IjawaObwb)- (8) 

mA.ms 


To account for our experimental conditions, we assume 
an isotropic distribution of 6 si /2 ground-state Cs atoms 
(nis = ± 1 / 2 ). (T+ polarized transitions thus lead to 

the excitation of four different states (f> corresponding to 
JA = JB = 3/2 and ttia = 1/2, 3/2 and tob = 1/2,3/2 
with equal weights. The intensity is obtained as the 
squared overlap p$( 0 , i?) = |(4>(i?)|(/)between the op¬ 
tically accessible function (p and the function 'I'(i?) ex¬ 
pressed in the molecule-fixed frame with the internuclear 
axis forming an angle 9 with the laboratory-fixed z axis. 
To calculate pqr{9, R), we transform the atomic functions 
jA.B^TiA.B into the molecule-fixed frame using the Wigner 
d rotation matrices and obtain 


p^{e,R) = \{^{RMf = 

I (^(^)l7AjAWA,7BjBWB) ^uJAmA (0)X 

jAmAjBmB 


(jATOA,jBTOB| 


( 9 ) 


We average p^(0, R) over all orientations 6 of the molec¬ 
ular axis relative to the z axis and obtain pi^(R) = 
/ p,s,{9, R)2TTsin{9)d6, which is directly proportional to 
the spectral intensity resulting from excitation of a 
molecular pair state jd/) at the interatomic separation 
R. The corresponding spectral line profile s^{E') is then 
determined as 


Sii'(£') — 

6 {E'- E^{R))j^-^R4^^M{R)R^dR, 

( 10 ) 
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FIG. 10 . (Color online) Calculated potential-energy curves for long-range Rydberg-atom pairs in the vicinity of np3/2np3/2 
asymptotes at a) n = 23 , b) n = 25 , and c) n = 32 and for = 0 , 1 , 2 , 3 (black, red, blue, and green curves, respectively) as 
obtained from our long-range potential model (details in the text). The intensity of the color reflects the value of p^{R), gray 
color denoting 0 % overlap with the atomic np3/2U.P3/2 state, full color denoting 1 % or higher overlap. The simulated spectra 
plotted in arb. units on the right-hand side of each of the three panels were calculated using Equation ( 111 . 


where Watom is the Rabi frequency of the atomic np 3/2 f— 
6 si /2 transition, Epp the asymptotic energy of the 
^P 3 / 2 tip 3/2 state, and the term 4a;ftom/(-E'pp “ E^s,{R))'^ 
accounts for the two-photon excitation of the pair- 
state with an off-resonant intermediate 'UP 3 / 26 S 1/2 (or 
6 si/ 2 np 3 / 2 ) state at a detuning of ^ (i?pp — E\i,{R)) p^ . 
The delta function ensures resonance with the molecu¬ 
lar state with the energy E^i,{R) and intensity p^{R). 
Assuming a gas of homogeneous density, the number of 
atom pairs at a separation R is proportional to 47 ri?^. 
For comparison with experimental spectra, the line pro¬ 
file s^(i?') is convoluted with the laser line profile G{E) 
and summed over all molecular states Idt) to obtain the 
full spectrum s{E) as 


/ OO 

s^{E') ■ G{E - E')dE' 

-00 


= IbTT^ 


^>00 /•OO 


6 {E' -E^{R)) 




' —00 J 0 


{Epp-E^(R)Y' 


p^{R)R^ ■ G{E - E')dRdE' 


= 167r^ 


0 {Epp — E^!{R)) 


.G{E - E^{R))x 


p^},{R)R^dR. 


( 11 ) 


The laser line profile G{E) is assumed to be a Gaussian 
function with the experimentally determined full width 
at half maximum of 140 MHz. The atomic Rabi frequen¬ 
cies Watom of the single-photon np 3/2 6 si /2 transitions 

act only as distance- and energy-independent overall scal¬ 


ing factors and are not taken explicitly into account when 
fitting the global amplitude of the simulated spectra s{E) 
to the experimental data (see Figs.an d§. 

The intensity of the color used to draw the potential- 
energy functions in Fig.[^ represents the value of the 
function |Hf(i?), gray corresponding to a negligible contri¬ 
bution of the optically accessible n.p 3 / 2 np 3/2 component 
and full color to a n.p 3 / 2 np 3/2 contribution of at least 1 %. 


At large distances, the only nonvanishing contributions 
are located in the immediate vicinity of the np 3 / 2 up 3/2 
asymptotes. At shorter distances, many more poten¬ 
tial curves contribute to the spectral intensity. How¬ 
ever, most potential curves are strongly attractive and 
the large gradients spread the overall intensity over 
broad spectral ranges, which renders the observation of 
pair states impossible. From the general appearance of 
the curves depicted in Fig.[T^ one expects sharp spec¬ 
tral structures to arise from regions where the potential 
curves that possess some np 3 / 2 np 3/2 character are flat. 
Such regions of enhanced intensity are encountered where 
the curves have a local maximum (a clear example of 
such a situation is given by the red and black curves cor¬ 
related to the 23si/224spi/2 asymptote in the left panel 
of Fig. 10 1 or closely approach their asymptotic value (all 
three curves correlated to nsi/ 2 (n + l)si /2 at large R val¬ 
ues in Fig. 10 1 . 

Sharp spectral structures are also encountered in re¬ 
gions where avoided crossings reduce the spectral den¬ 
sity, leading to intensity minima. Characteristic exam¬ 
ples of this situation are the minima observed in the 
high-frequency wings of the n-P 3 / 2 /iP 3/2 resonances for 
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FIG. 11. (Color online) Calculated potential-energy func¬ 
tions near the 24p3/2 24p3/2 pair-dissociation asymptote (left 
panel), simulated spectrum s{E) (middle, arb. units), and 
experimentally measured Cs''"-ion signal as a function of the 
two-photon detuning 2Ai> (right panel, arb. units). The color 
code is the same as in Fig. |10| The horizontal dashed lines 
indicate the positions of minima and maxima in the simulated 
spectrum for comparison with the experimental spectrum and 
the calculated potential-energy functions. 


n = 22 — 26 in Fig.[^ 


B. Interpretation of the UV-laser spectra 


The procedure we followed to interpret the exper¬ 
imental UV-laser spectra of Rydberg-atom pairs such 
as those presented in Figs. BM and was to evaluate 
the spectrum s{E) using Eq. (11) retaining successively 
higher terms of the long-range interaction series (Eq. if), 
which, in turn, necessitated the extension of the basis set 
I’T'A^aJawA) JbWb) to include pair states with Za and 
Zb values up to 6. Convergence of the eigenvalue calcu¬ 
lations was checked by increasing the range of principal 
quantum numbers riA and ub- Because the atomic ra¬ 
dial transition matrix elements rapidly decrease with the 
value of An, only states differing in effective principal 
quantum number n* by no more than 1 were considered 
in the calculations. It further turned out to be suffi¬ 
cient to only consider the pair states with asymptotes 
located within an energy window of 0.64 • 4:Rcs/{n*)^ 
around the relevant np 3 / 2 ?T'P 3/2 asymptote. States lo¬ 
cated outside this range only approach the np 3 / 2 ?T-P 3/2 
potential curves at short range, where the gradients of 


the potential-energy curves are large and spectral struc¬ 
tures are broad and weak. 

As explained in Section IV A, the model does not in¬ 
clude adjustable parameters, but relies on the known 
quantum defects of the Cs-atom Rydberg states. The 
structures of the spectra in the regions located below the 
np 3 / 2 np 3/2 dissociation asymptote (see Fig.j^ are pri¬ 
marily determined by the dipole-dipole interaction term 
which couples the np 3 / 2 np 3/2 and upi/ 2 ?^P 3/2 pair states 
with the nsi/ 2 (n + l)si /2 pair states. The potential- 
energy functions that give the nSi/ 2 (ri+l)si /2 resonances 
their characteristic interaction-induced shifts (see inset 
of Fig.|4| and fine structures (see Fig.jsl) are depicted in 
Fig.[T0' for n = 23, 25 and 32. The calculated spectra 
s{E) are depicted on the right-hand side of each panel, 
and dashed horizontal lines highlight the origins of the 
sharp spectral structures. 

Three states correlating to the nsi/ 2 (n-|-l)si /2 asymp¬ 
totes possess a significant admixture of the optically ac¬ 
cessible np 3 / 2 rrp 3/2 pair state, two of which are U = 0 
(black lines) and the third (red line) is a U = 1 state. 
At the highest n values studied experimentally (27 < 
n < 37), for which the results obtained at n = 32 are 
representative (see Fig. 101, all three states are attrac¬ 
tive. In the range of internuclear distances where the 
np 3 / 2 np 3/2 character is significant, the dipole-dipole in¬ 
teraction with the n.p 3 / 2 ?T'P 3/2 state is the dominant in¬ 
teraction for all three states. The interaction with the 
^Pi/ 2 ’^P 3/2 pair state, which is located more than 10 GHz 
below the np 3 / 2 »^P 3/2 asymptote at n = 32, does not 
significantly affect the nsi/ 2 (n + l)si /2 potential curves 
in this range of n values. The calculations predict a 
red detuning of the pair resonance with respect to the 
nSi/ 2 (rr+l)si /2 asymptote and a red-degraded line shape, 
in excellent agreement with the experimental observa¬ 
tions (see Fig.|^. 

As n decreases below n = 28, the position of the 
nsi/ 2 (n + l)si /2 asymptotes rapidly shifts away from 
the 'np 3 / 2 np 3/2 asymptote and toward the npi/ 2 np 3/2 
asymptote. Consequently, the effect of the dipole-dipole 
interaction with the ?rpi/ 2 *rp 3/2 pair state becomes im¬ 
portant: The H = 1 and one of the two H = 0 potential 
functions become repulsive at intermediate distances be¬ 
yond about 0.3 /xm at n = 25 and beyond about 0.20 /xm 
at n = 22. At shorter range, the interaction with the 
rrP 3 / 2 riP 3/2 pair state is always dominant and makes the 
nsi/ 2 (n + l)si /2 potentials attractive, resulting in local 
maxima in these two curves. These two maxima give 
rise to two peaks in the calculated spectra located on 
the high-frequency side of the nsi/ 2 {'n + l)si /2 asymp¬ 
totes. The third peak, which corresponds to the position 
of the asymptote, is well resolved below n < 25. The 
calculated spectra s{E) accurately reproduce the exper¬ 
imental spectra at all n values investigated (see Fig.|^. 
The interaction-induced shifts of the pair resonances ex¬ 
tracted from simulated and experimental spectra agree 
within the experimental uncertainty, as shown in the 
inset of Fig. Even the positive interaction-induced 
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shift of the n = 25 pair resonance with respect to the 
25 si/226si/ 2 asymptote (see Fig.® is accounted for by 
the calculations, which indicate tnat it results from the 
overlap of signals originating from the long-range part of 
the potentials and from the flat potential region near the 
maximum of the 17 = 1 curve (red trace in Fig. |10[ 3)). 

The potential-energy functions of the 17 = 0 (black), 
1 (red), 2 (blue), and 3 (green) states required to 
interpret the spectral structures located on the high- 
frequency side of the 24p3/224p3/2 asymptotes are de¬ 
picted in Fig. El Many pair states correlating to the 
20f5/2,7/225Si/2, 20f5/2 J/221f5/2,7/2, 22d3/224d3/2,5/2j 

21gj24si/2, 20gj25si/2 and 21hj24si/2 asymptotes are lo¬ 
cated in this region, but they only possess 24p3/224p3/2 
character near avoided crossings with the np 3 / 2 ?^P 3/2 
(17 = 0 — 3) states. In the absence of these crossings, 
the repulsive curves associated with the 24p3/224p3/2 
asymptote would lead to a broad and structureless high- 
frequency wing of the atomic transition. At the avoided 
crossings no curve possesses significant 24p3/224p3/2 
character. Consequently, the spectral excitation func¬ 
tion s{E) has minima that can be directly linked to the 
avoided crossings, as indicated by the horizontal dashed 


lines in Fig. 11 


The simulations of the spectra presented in Fig. [^indi¬ 
cate that a potential model restricted to the dipole-dipole 
interaction fails to reproduce the details of the spectral 
structures (left panel) on the high-frequency side of the 
atomic transitions and that the successive inclusion of 
interaction terms up to the octupole-octupole term grad¬ 
ually improves the agreement between measured and cal¬ 
culated spectra. 

The fact that the calculations almost perfectly repro¬ 
duce the experimental spectra (see Figs.[^ and is re¬ 
markable and clearly reveals that (i) the long-range adi¬ 
abatic potential model is adequate both with respect to 
the multipole expansion and the basis set, and (ii) pair¬ 
wise interactions are sufficient to fully describe the be¬ 
havior of the Rydberg gas at the time scale probed by 
our laser pulses (« 5 ns) under the experimental condi¬ 
tions used to record the spectra. 

It may appear surprising at hrst sight that adiabatic 
potential-energy functions can be used to accurately de¬ 
scribe the spectra of Rydberg-atom pairs. Indeed, the 
high density of electronic states implies slow electronic 
dynamics, which in general tends to be accompanied by 
a breakdown of the Born-Oppenheimer and adiabatic 
approximations. The reason why the potential-energy 
curves provide an accurate description of the spectra is 
that the electronic dynamics, though slow, is still several 
orders of magnitude faster than the nuclear dynamics. 
For any given value of 17, the state density of electronic 
Rydberg-atom pairs is approximately 1 — 5GHz“^ at 
u = 24 which implies that the timescale of the electronic 
motion is about 1 ns. At the translational temperature of 
40 /rK of the Cs vapor, the mean velocity of the Cs atoms 
is about 10 cm/s, so that they can be considered station¬ 
ary on the timescale of the electronic motion, and also 


during the UV-laser pulse. At internuclear distances be¬ 
yond R = 200 nm, the gradients of the potential curves 
are less than 500 GHz/p,m so that the relative velocity 
acquired by the atoms in Ins is less than 2m/s which 
is also negligible. Over longer timescales, the accelera¬ 
tion of the atoms in the attractive potentials inevitably 
leads to nonadiabatic effects and eventually to Penning 
ionization. 


C. Simulation of Penning ionization 

To analyze the Penning-ionization experiments de¬ 
scribed in Section HID (see in particular Fig.^), we use 
a simple model which relies on the spectral excitation 
function s{E) and treats the relative motion of the atoms 
on the two 17 = 0 and the 17 = 1 potential curves asso¬ 
ciated with the 32 si/233si/ 2 asymptote classically. All 
other pair states are excited at short range (i.e., at R 
values below 0.4/rm) and only contribute to the spec¬ 
tra weakly on the low-frequency side of the resonance. 
The potential curves of the relevant molecular states are 
drawn in color in Fig.j^). Below 0.7/im, the potential- 
energy curves reveal several avoided crossings. Following 
photoexcitation, the two interacting Rydberg atoms are 
accelerated towards each other until Penning ionization 
takes place, which is assumed, for simplicity, to happen 
at an interatomic distance of 200 nm. The probabilities 
of adiabatic traversals of the avoided crossings are esti¬ 
mated with the Landau-Zener formula 

Tdiab = exp | —2-77I ; F’adiab = 1 — Fhiab, (12) 

where V 12 represents the minimal distance between the 
potential curves and the slew rate dE/dt is determined 
from the relative motion of the two atoms on the po¬ 
tential curves. The first avoided crossing of the 17 = I 
state (see red curve in Fig.[^)) is traversed adiabatically 
(Fadiab > 99%). The adiabatic potential has a minimum 
at i? « 0.58/im and the state therefore does not rapidly 
reach the region where Penning ionization takes place. 
The motion on the two 17 = 0 potential curves (black 
and blue/green lines in Fig.[^)) is such that the regions 
of short internuclear distances where Penning ionization 
takes place is reached rapidly despite several avoided 
crossings. The simulation of the ionization dynamics of 
the 17 = 0 states adequately reproduces both the rise of 
the spontaneous-ionization signal and the gradual shift of 
the peak position (see insets of Fig.[^)) observed exper¬ 
imentally for delay times up to « 2 /is. For longer delay 
times a slower rise of the spontaneous-ionization signal 
up to « 10 /ts is observed, which is not explained by our 
model (see lower inset of Fig.j^)). We attribute this ad¬ 
ditional ion signal to the ionization of pairs of atoms in 
the 17 = 1 state. A possible mechanism could be a weak 
mixing of 17 = 0 and 17 = 1 states {e.g. by a resid¬ 
ual electric field), which would lead to a small avoided 
crossing of the two symmetries (see red and green/blue 
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traces in the inset of Fig. 8b) at i? ~ 0.5 /rm) and open 
a path for ionization of pairs in the metastable 14 = 1 
state. The probability for an adiabatic crossing is in¬ 
versely proportional to the velocity of the two colliding 
atoms at the crossing (see Eq. For pairs excited at 

lower frequencies, this velocity decreases and the ioniza¬ 
tion probability thus increases. This is consistent with 
the observed increase of the deviation between experi¬ 
mental and simulated ionization signals at longer delays 
when the excitation laser is tuned to lower frequencies 
(see lower inset in Fig. 8a)). For long delay times, the 
experimental spontaneous-ion signal approaches approxi¬ 
mately 50% of the initial total-ion signal, which indicates 
that every pair of atoms eventually undergoes a Penning 
ionization process. 

The model discussed in this section predicts the cor¬ 
rect qualitative behavior of fast and slow spontaneous 
ionization processes, but does not satisfactorily explain 
the observation of an almost constant total-ion signal for 
all experimental delay times (see dashed-colored traces in 
Fig. [8^)). Indeed, we predict that the Rydberg atom un¬ 
dergoing a transition to a lower Rydberg state in the Pen¬ 
ning ionization process must end in a level with n* <21, 
whereas the amplitude of the pulsed electric field used for 
field ionization (1250 V/cm) should only ionize Rydberg 
states with n* > 22. 


D. Simulation of millimeter-wave spectra between 
Rydberg-atom—pair states 

The potential model and the spectral distribution func¬ 
tion s{E) presented in Section IV A were also used to cal¬ 
culate the millimeter-wave spectra of the 32si/233p3/2 t— 
3283 ^ 2 ^ 13 ^ 1/2 Sjud 3283 ^ 2 ^ 383^2 —^ 31p3^233s3y2 transi¬ 
tions displayed in Figs. [^) and c). To this end, we 
used, as initial distribution of Rydberg-atom-pair states 
the relative population in the two 12 = 0 and the 
12 = 1 states that are accessed optically at the peak 
of the 32 s3/233s 3/2 resonance, and consider transitions 
from these states to the pair states correlating to the 
31p3/233s3/2 and 32s3/233p3/2 asymptotes that are ac¬ 
cessible in a single-photon millimeter-wave transition. 
The relevant potential-energy functions are displayed in 

Fig.[Tg 

The attractive nature of the three potential func¬ 
tions associated with the 32s3/233s3/2 asymptote im¬ 
plies, in a classical description of the nuclear mo¬ 
tion, that the internuclear separation decreases with 
time. Because the potential curves of the states asso¬ 
ciated with the 31p3/233s3/2 and 32s3/233p3/2 asymp¬ 
totes are less attractive than those of the states asso¬ 
ciated with the 32 s3/233s 3/2 asymptote, the frequency 
of the 32 s3/233s 3/2 —t 31p3/233s3/2 transition decreases, 
and that of the 32s3/233p3/2 ^ 32s3/233s3/2 increases, 
with time. This evolution (chirp) of the transition fre¬ 
quencies reduces the transition probabilities. The narrow 
bandwidth (1 MHz) of the millimeter-wave radiation im- 



N 
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FIG. 12. (Color online) Potential-energy curves of the 
states correlated to the 32si/2 33si/2 (middle panel), the 
32si/233p3/2 (top panel) and the 31p3/233si/2 (bottom panel) 
asymptotes used for the simulation of the millimeter-wave 
spectra presented as red dashed lines in Figs. and c). 
Black (grey), red and blue designate states with f2 = 0,1 and 
2, respectively. The intermediate panels present the transi¬ 
tion probabilities calcnlated for the absorption (Ptrans t) a-nd 
stimulated emission (Ptrana 4) out of the 12 = 0 (black) and 
12 = 1 (red) states correlated to the 32si/233si/2 asymptote. 
All energies are given relative to the energy of the 32si/2 33si/2 
asymptote. See text for details. 


plies that the main contribution to the millimeter-wave 
spectra originate from pairs with large internuclear sepa¬ 
ration R, for which the evolution (chirp) of the transition 
frequency is slow. 

To model the intensity distribution, we introduce a 
cut-off chirp of 2 MHz/ps, beyond which we assume that 
the transitions are too weak and broad to be observed 
and determine the range of internuclear separations con¬ 
tributing to the millimeter-wave transition intensities ac¬ 
cordingly. In this range, we calculate the distribution of 
transition frequencies from the potential-energy curves. 
For each frequency, we determine the spectral intensities 
as sum over products of the initial populations in the 
12 = 0 and 1 states correlated to the 32s3/233s3/2 asymp¬ 
tote and the relative transition probabilities 

Ttrans(R) = 

( 13 ) 
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where is the cutoff function described above, 

Ao,±i is the sum of cr^ and tt electric-dipole transition op¬ 
erators, and the double sum runs over all initially excited 
pair-state d'i for a fixed value of and all dipole-coupled 
final states 'h/. 

The calculated spectra, convoluted with the experi¬ 
mental line-shape function of the corresponding atomic 
transitions, are displayed as dashed red lines in Figs. [6>) 
and c) and are in qualitative agreement with the experi¬ 
mental ones. In particular, the shifts from the asymp¬ 
totic transition frequencies, which correspond to the 
spectra of the isolated atoms displayed as lower traces, 
and the asymmetric line shapes (red-degraded for the 
32 si/233si/2 —>■ 31p3/233si/2 transition and blue de¬ 
graded for the 32si/233p3/2 ^ 32si/233si/2) are well 
described by the calculations. This agreement repre¬ 
sents a further illustration of the ability of our potential 
model to describe the spectral and dynamical properties 
of Rydberg-atom pairs. 


V. CONCLUSIONS 

In this article we have presented a comprehensive sur¬ 
vey of the Rydberg excitation spectrum in a dense (den¬ 
sity 2 • 10^^ atoms/cm^), ultracold (T = 40 /rK) Cs gas 
recorded using an intense (peak intensity 10® W/cm^) 
narrow-band (bandwidth 140 MHz) UV laser source. 

Next to the strong lines corresponding to the atomic 
^P 3/2 6s transitions, the spectra also exhibit a very 
rich structure of weaker, sharp spectral features, which 
are attributed to pairs of Rydberg atoms, also known as 
Rydberg macrodimers [101117], interacting through long- 
range electrostatic forces. Sharp features of three differ¬ 
ent origins have been identified: (i) lines with maxima 
at positions located very close to nsi/ 2 (n -I- l)si /2 and 
nsi/ 2 (n — 3)f dissociation asymptotes and which corre¬ 
spond to the formation of weakly interacting Rydberg- 
atom pairs at long range, (ii) lines with maxima located 
above the nsi/ 2 (n -I- l)si /2 dissociation asymptotes and 
which correspond to excitation to flat potential regions 
near local maxima of the potentials of H = 0 and 1 
molecular states correlated with nsi/ 2 (n -I- l)si /2 disso¬ 
ciation asymptotes, and (iii) ’’window” resonances lo¬ 
cated on the high-frequency side of the atomic transi¬ 
tions, which arise from avoided crossings between repul¬ 
sive molecular potential curves correlated to np 3 / 2 np 3/2 
dissociation asymptotes and attractive potentials near 
(n — 4)fj(n — 3)fj/ (j, j' = 5/2, 7/2) dissociation asymp¬ 
totes. The first category of resonances had been observed 
in earlier studies (see, e.g. Refs. HTHIll EUEg). The 
observation of the latter two categories (ii) and (iii) rep¬ 
resents the first major result of our investigation. 

The second major result is the demonstration, on the 
basis of a rigorous comparison of experimental and calcu¬ 
lated spectra, that long-range molecular potential models 
provide a quantitatively correct and accurate description 
of the behavior of interacting Rydberg atoms in ultra¬ 


cold gases, provided that a sufficient number of channels 
and long-range interaction terms are considered. The po¬ 
tential model derived to describe Cs Rydberg-atom pairs 
adequately reproduces experimental UV-laser photoas¬ 
sociation spectra, millimeter-wave spectra of transitions 
between different Rydberg-atom-pair states, and the ob¬ 
served Penning-ionization dynamics of the Rydberg-atom 
pairs. Given the very high density of electronic states, it 
first appears astonishing that adiabatic potential models 
are at all adequate to describe the spectral and dynam¬ 
ical properties of interacting Rydberg-atom pairs. The 
reason for the usefulness of such models is that, in the 
range of interaction strengths and internuclear distances 
that are relevant for most experiments in ultracold Ryd¬ 
berg gases, the heavy nuclei are almost stationary on the 
time scale associated with the electronic motion, taken 
as the inverse state density, i.e., about 1 ns at n « 24. 

In the discussion of dipole- and van der Waals-blockade 
effects in the excitation of Rydberg states, the attractive 
and repulsive potentials arising from the dipole-dipole 
interaction near the asymptotes of Rydberg-atom-pair 
states, are often represented, for clarity and simplicity, 
as two-level systems of the kind depicted in Fig. j^. The 
present analysis of Rydberg-atom-pair states of Cs in 
the range of principal quantum number n = 22 — 36, 
which includes all relevant potential curves and long- 
range interaction terms, demonstrates the limits of two- 
level Rydberg-blockade models. Indeed, many inter¬ 
acting molecular potentials give rise to a multitude of 
avoided crossings, resonances and level shifts. The range 
of observable phenomena is thus much wider than two- 
level blockade models predict, which on the one hand 
reduces their range of applicability, but on the other also 
offers new opportunities for scientific investigations and 
applications. 
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Appendix: Calculation of matrix elements of VL 

In this appendix, explicit expressions for the ma¬ 
trix elements of the atomic multipoles contributing to 
kinter (seC Eq. @) are derived in the product basis 
I^^A^AJAWA, 'raB^BjBWB). 

The matrix elements of an atomic multipole 

Ql,. (r-) = (r) (A. 1) 
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in this basis are given by 

{'n%jW^ QlA"^ \n^kj^uJ^) 


where are the Clebsch-Gordan coefficients and 


AtT t '' 

21 I ^ \n^kj^) YL^ujif) \kj,uJi ), 

(A.2) 

where the radial matrix elements 

\ndiji) = RZAai 

= J K'j'j'ir)r^'^nd,jAr)dr (A.3) 

are evaluated numerically using the Numerov method [45] 
and a model core potential |46| . The angular matrix ele¬ 
ments are evaluated using the analytical expressions [15] 

YlA^) 

= (^_iy.+i’+L+h 

" {U t} ’ <*■“> 


Jl J 2 J 3 

ji js je 


is Wigner’s 6-j symbol. 


As an example we give the full expression for the ma¬ 
trix elements of the lowest-order term in the multipole 
expansion, the dipole-dipole interaction: 


I^A^AjAd-’Ai ^B^BjBd.’B) 

_ Wa-I-Ib+^a+^b 7?“3 p^-bAIb 

- y ^ ^tiaIaJa ^nslujB 

^ \/(2ja -|- 1)(2jb + 1)(2/a + 1)(2?b + 1) 

w f ^Iai'^a ^ 1 b>“b I 

I o/^Iai'^a ^Ib’^b I ^JA’^a ^3b’‘^b a 
^'~'jA‘^A,10'~^jBi^B,W '-^JAl->A,l-l'^iB<.-JB,ll J 

P^aO f^A 2 f^B 2 /a r\ 
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